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Abstract 

First-principles  density  functional  theory  (DFT) 
calculations  using  the  generalized  gradient  approximation 
(GGA)  have  been  used  to  expand  the  analysis  of  the 
adsorption  properties  of  polynitrogen  and  high  nitrogen 
compounds  on  Al(lll)  and  AlN(OOOl)  surfaces.  The 
electron-ion  interaction  has  been  described  using  ultrasoft 
(PAW)  pseudopotentials  in  the  case  of  Al  (AIN).  All 
calculations  have  been  done  using  the  PW91  generalized 
gradient  approximation.  In  the  case  of  Al,  the  calculations 
employ  periodic  slab  models  with  4  layers,  ranging  in  size 
from  (4x4)  to  (6x6)  surface  unit  cells,  and  containing  up  to 
144  Al  atoms.  For  AIN,  slab  models  with  (2x2)  and  (3x3) 
unit  cells  and  with  8  an  4dual  AIN  layers,  respectively,  have 
been  considered.  Complementary  quantum  chemical 
calculations,  utilizing  DFT  and  second-order  perturbation 
theory  methods,  of  the  ground  state  potential  energy 
surfaces  of  the  corresponding  polynitrogen/high  nitrogen 
species  in  the  absence  of  the  aluminum  surface  also  have 
been  performed.  Previous  results  obtained  in  the  first  two 
years  of  this  challenge  project,  related  to  adsorption 
properties  of  all  nitrogen  Nx(x=l-6,8-12)  compounds  or  of 
non-substituted  and  mono-substituted  triazenes  have  been 
extended  to  the  case  of  high  nitrogen  heterocycle  systems. 
Specifically,  we  have  considered  representative  six- 
membered  P-N-C  heterocycles  containing  azide  ligands  on 
the  C  and  P  atoms  such  as  C2N15P,  CNigPz  and  N21P3 
systems.  Additionally,  we  extended  our  analysis  by 
including  compounds  containing  the  nitrogen  rich  CNf 
anion  (i.e.,  the  energetic  salts  NHfiCNf,  N2H5+CN7',  and 
CN4H7+CN7~)  and  heterocyclic-based  ionic  systems 
(C2N3H4+-  NO 3 )  for  which  direct  optimizations  and  ab 
initio  molecular  dynamics  simulations  have  been 
performed.  Finally,  the  interaction  properties  of  atomic 
and  molecular  oxygen  and  of  other  energetic  species  like 
nitromethane  ( CH3NO? )  with  AlN(OOOl)  and  A/iV(0001) 
surfaces  have  been  analyzed.  For  these  systems  the 
individual  bonding  mechanism  involving  either 
nondissociative  ( molecular )  or  dissociative  processes  has 
been  determined,  depending  on  both  the  molecular 
orientation  and  the  surface  sites  involved.  The  binding 
energies  were  found  to  be  highly  sensitive  to  the  specific 
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adsorbed  species  and  its  local  configuration.  Dissociation 
processes  were  also  found  to  be  highly  exothermic. 

1.  Introduction 

Investigation  of  the  interactions  of  energetic  high- 
nitrogen  compounds  with  ultrafine  or  nanophase  aluminum 
particles  is  a  topic  of  current  DoD  interest  with  potential 
applications  in  rocket  and  missile  propulsion.  One  of  the 
topics  of  prime  importance  is  to  determine  if  the  thermal 
instability  inherent  in  many  high  nitrogen  compounds  might 
be  reduced  by  adsorption  or  chemisorption  onto  the  surface 
of  aluminum,  while  simultaneously  forming  a  protective 
coating  on  the  metal  surface  which  minimizes  the  formation 
of  an  inert  oxide  surface  layer. 

2.  Computational  Methods 

The  calculations  performed  in  this  study  were  done 
using  the  Vienna  ab  initio  simulation  package  (V  ASP). 11-31 
This  program  evaluates  the  total  energy  of  periodically 
repeating  geometries  based  on  density-functional  theory  and 
the  pseudopotential  approximation.  For  calculations  on  the 
Al  surface  the  electron-ion  interaction  is  described  by  fully 
non-local  optimized  ultrasoft  pseudopotentials141  while  the 
projector  augmented  wave  (PAW)  method  of  Blochl  151  was 
used  for  the  AIN  surface.  This  difference  in  treatment  was 
due  to  availability  of  PAW  pseudopotentials  for  partially 
charged  H  atoms  used  in  the  AIN  slab  calculations.  Periodic 
boundary  conditions  were  used,  with  the  one -electron 
pseudo-orbitals  expanded  over  a  plane-wave  basis  set  with 
cutoff  energies  of  435  eV  for  Al  and  500  eV  for  AIN.  These 
cutoff  energies  are  consistent  to  the  precision  level  required 
by  the  VASP  code  for  the  chemical  species  of  interest. 
Calculations  were  performed  using  the  generalized  gradient 
approximation  (GGA)  density  functional  theory  with  the 
PW91  exchange-correlation  functional.161  The  sampling  of 
the  Brillouin  zone  was  performed  using  a  Monkhorst-Pack 
scheme.171  The  minimum  energy  path  between  different 
minima  was  determined  using  the  climbing  image  nudged 
elastic  band  method  developed  by  Jonsson  and  co¬ 
workers.18'91 


This  page  is  Distribution  A:  approved  for  public  release;  distribution  unlimited. 


The  GAMESS1 10  quantum  chemistry  code  was  used  to 
predict  the  structures  and  properties  of  isolated  polynitrogen 
and  high-nitrogen  compounds.  The  theoretical  methods 
employed  included  second  order  perturbation  theory1111 
(MP2,  also  known  as  MBPT(2))  and  density  functional 
theory  (DFT)  using  the  hybrid  B3LYP1121  functional.  The 
6-31  lG(2df,p)  and  6-311++G(2df,p)1131  basis  sets  were  used 
for  nonionic  and  ionic  species,  respectively. 

3.  Results  and  Discussion 

Extension  of  our  comprehensive  analysis  of  the 
interaction  between  high-nitrogen  and  polynitrogen 
compounds  with  aluminum  surface  was  done  by 
considering  four  new  classes  of  nitrogen  systems. 
The  first  class  consists  of  high  nitrogen  heterocyclic 
systems.  In  this  case  we  have  considered  the  case  of  six- 
membered  P-N-C  heterocycles  containing  azide  ligands  on 
the  C  and  P  atoms  such  as  C2Ni5P,  CNi8P2  and  N21P3. 

The  second  class  is  comprised  of  salts  of  the  CN7’  anion. 
As  shown  very  recently,1141  such  systems  represent 
promising  candidates  for  development  of  nitrogen-rich 
green  energetic  materials.  For  this  set  of  compounds  we 
have  analyzed  several  high  explosive  salts  containing  the 
CN7"  anion,  namely  the  hydrazinium  ([N2H5][CN7]), 
ammonium  ([NH4][CN7])  and  aminoguanidinium 
([CN4H7][CN7])  salts.1141 

The  third  class  of  systems  includes  compounds  which 
combine  an  oxygen-rich  dinitramide  anion  with  a  nitrogen- 
rich  cation  such  as  l,5-diamino-4-methyltetrazolium 
dinitramide  ([C2N6H6][N304]).  Such  systems1151  represent 
promising  candidates  for  development  of  new  tetrazole- 
based  oxidizers.  Finally,  we  extended  our  analysis  to  the 
case  of  heterocyclic -based  ionic  systems  (C2N3H4+-N03’) 
for  which  direct  optimizations  and  ab  initio  molecular 
dynamics  simulations  have  been  performed. 

The  adsorption  studies  of  the  above  described  systems 
have  been  studied  using  3D  periodic  slab  models  of 
different  sizes.  These  range  from  (4x4)  surface  supercells 
used  in  the  case  of  small  systems  to  (6x6)  surface  units  for 
the  largest  systems.  In  all  cases,  the  slab  models  contained 
four  layers  of  atoms.  The  optimizations  have  been 
performed  for  both  the  adsorbate  and  the  A1  atoms  in  the 
top  two  layers  while  the  bottom  two  layers  were  frozen  at 
the  bulk  optimized  conditions.  Adsorption  at  different 
surface  sites  and  for  different  molecular  orientations  have 
been  investigated  and  the  corresponding  adsorption  energies 
have  been  calculated  using  the  expression 
Eads- Enioiec"i"Esiafo— E(moiec+siab>,  where  Em0]ec  is  the  energy  of 
the  isolated  molecular  species  at  its  equilibrium  geometry, 
Eslab  is  the  total  energy  of  the  isolated  aluminum  slab,  and 
E(moiec+siab)  is  the  total  energy  of  the  adsorbate/slab  system. 
A  positive  Eads  corresponds  to  a  thermodynamically  stable 
adsorbate/slab  system. 

P-N-C  Heterocycles  Adsorption  on  on  AI(lll).  The 

adsorption  properties  of  C2N15P,  CN18P2  and  N2iP3  systems 


have  been  evaluated  using  a  (6x6)  slab  model. 
Representative  adsorption  configurations  are  presented  in 
Figure  1  and  the  corresponding  binding  energies  are 
summarized  in  Table  I.  In  the  case  of  C2N15P,  several 
distinct  binding  configurations  were  found,  including  non- 
dissociative  adsorptions.  A  representative  configuration  of 
this  type  is  shown  in  Figure  la  in  which  adsorption  takes 
place  with  formation  of  multiple  Al-N  bonds.  Additional 
similar  configurations  were  found,  with  adsorption  energies 
ranging  from  60.1  to  65.9  kcal/mol.  The  presence  of  N3 
groups  in  close  proximity  to  the  surface  can  also  lead  to 
bond  fragmentation  and  elimination  of  one  (Figure  lb)  or 
two  (Figures  lc  and  Id)  molecules  of  N2.  The  dissociative 
adsorption  of  C2N15P  was  found  to  be  highly  exothermic, 
with  dissociation  enthalpies  as  high  as  92.3  kcal/mol 
relative  to  the  adsorbed,  undissociated  state  (Figure  la). 

Similar  adsorption  characteristics  were  found  for  the 
other  two  heterocycle  systems  CNi8P2  and  N2iP3. 
Specifically,  molecular  (undissociated)  adsorption  states 
were  found  to  be  possible  for  both  types  of  species  (see 
Figures  le,  li,  and  lj).  The  non-dissociative  configurations 
form  varying  number  of  Al-N  bonds,  which  in  turn  lead  to  a 
spread  in  the  calculated  binding  energies.  For  example,  the 
binding  energies  of  the  non-dissociative  configurations 
illustrated  in  Figures  li  and  lj  are  49.5  kcal/mol  and  80.9 
kcal/mol,  respectively,  in  which  four  and  six  Al-N  bonds, 
respectively,  are  formed.  Dissociative  adsorption  was  also 
found  to  take  place  for  both  CNi8P2  and  N21P3  compounds 
and  the  overall  processes  are  highly  exothermic. 
Representative  dissociative  binding  configurations  are 
illustrated  in  Figure  1,  corresponding  to  formation  of 
CN16P+N2,  (panel  If),  CN14P  +2N2  (panels  lg  and  lh), 
N19P3+N2  (panel  lk),  and  N17P3+2N2  (panel  11).  In  all  these 
cases  formation  and  desorption  of  N2  molecules  was  found 
to  be  thermodynamically  favored.  The  presence  of  several 
Al-N  bonds  in  these  cases  contributes  to  the  increased 
stability  of  these  systems. 

Adsorption  of  CN7  Salts  Systems.  A  primary 
objective  of  the  adsorption  studies  of  [CN7]"  salts  on  the 
Al(l  1 1)  surface  is  to  determine  how  ionic  adsorption  differs 
from  adsorption  of  neutral  compounds.  For  all  three  cases 
investigated;  i.e.  for  hydrazinium  ([N2H5][CN7]), 
ammonium  ([NH4][CN7])  and  aminoguanidinium 
([CN4H7][CN7])  salts,  the  anion  was  found  to  bind  directly 
to  the  A1  surface.  Additionally,  the  anion  serves  as  an 
anchor  center  for  the  cation  system.  Such  representative 
adsorption  configurations  are  depicted  in  Figure  2. 

In  the  case  of  hydrazinium  5-azidotetrazolate,  the 
[CN7]"  anion  can  adsorb  either  non-dissociatively  (see 
Figure  2a)  or  dissociatively  with  elimination  of  a  N2 
molecule  (see  Figure  2b).  In  both  instances,  bonding  to  the 
surface  is  mediated  by  formation  of  several  Al-N  bonds. 
The  resultant  CN7  or  CN5  fragments  bonded  directly  to  the 
surface  also  function  to  anchor  the  corresponding  [N2H5]+ 
cation.  Furthermore,  this  can  take  place  either  with  (Figure 
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2c  and  2d)  or  without  (Figures  2a  and  2b)  proton  transfer 
from  the  cation  to  the  anion,  with  the  former  process  being 
energetically  favored  (see  Table  1.)  Additionally, 
adsorption  of  the  [CN7]“  anion  with  subsequent  elimination 
of  N2  molecules  is  also  more  exothermic  than  non- 
dissociative  adsorption. 

In  the  case  of  ammonium  5-azidotetrazolate 
([NH4]+[CN7]";  see  panels  2e-2h)  the  primary  bonding  to  the 
surface  is  also  through  the  [CN7]  anion.  Not  surprisingly, 
the  corresponding  adsorption  energies  increase  with  the 
number  of  Al-N  bonds.  The  cation  remains  bonded  to  the 
adsorbed  anion  and  in  all  instances  charge  neutralization 
takes  place  via  proton  transfer  from  the  cation  to  the  anion. 
The  energetically  most  favorable  process  is  dissociative 
absorption  of  the  anion,  leading  to  desorption  of  N2  and 
chemisorption  of  the  remaining  CN5  fragment,  as  illustrated 
in  Figure  2h. 

Similar  behavior  is  observed  in  the  case  of  the  1- 
aminoguanidinium  5-azidotetrazolate  ([CN4H7]+[CN7]")  salt 
(see  Figures  2i-21.)  The  adsorption  energies  increase  with 
the  increase  in  number  of  Al-N  bonds  between  [CN7]  and 
the  A1  surface.  Additionally,  dissociative  adsorption  of  the 
anion  with  formation  of  N2  is  found  to  be  energetically 
favorable.  A  unique  feature  of  this  system  is  that  no 
proton  transfer  from  the  [CN4H7]+  cation  to  the  [CN7]" 
anion  is  observed  for  any  of  the  adsorption  configurations 
investigated. 

Adsorption  of  l,5-diamino-4-methvltetrazolium 
dinitramide  (DMTD)  In  the  case  of  DMTD,  which 
contains  both  an  oxygen-rich  dinitramide  anion  and  a 
nitrogen-rich  cation,  adsorption  to  the  surface  can  take 
place  through  both  the  anion  and  the  cation.  In  the  majority 
of  cases  investigated,  adsorption  takes  place  with 
dissociation  of  dinitramide  anion  (see  Figure  3)  and 
formation  of  strong  Al-O  bonds  via  isolated  O  atoms  or 
N02  fragments.  Additionally,  formation  and  desorption  of 
N2  molecules  is  also  observed.  The  overall  reactions  are 
highly  exothermic  due  to  the  large  energy  release  upon 
formation  of  Al-O  bonds. 

Adsorption  of  O  species  on  AIN  Surface  One  of  the 

fundamental  issues  considered  in  this  study  is  the  possibility 
of  using  AIN  as  a  coating  to  inhibit  or  prevent  formation  of 
oxide  layers  on  the  aluminum  surface.  For  this  purpose  the 
adsorption  properties  of  atomic  and  molecular  oxygen  on 
different  surfaces  of  AIN  are  analyzed.  Given  the  nature  of 
the  AIN  crystal  with  alternating  layers  of  A1  and  N  atoms, 
both  the  A1  terminated  AlN(OOOl)  surface  and  the  N 
terminated  AIN  (0001 )  surface  are  studied.  These 
calculations  have  been  done  using  2x2  slab  models  with 
eight  AIN  layers  and  terminated  with  fractionally  charged  H 
atoms.  The  role  of  these  pseudohydrogens  is  to  saturate  the 
dangling  bonds  in  the  bottom  of  the  slabs  and  to  prevent 
unphysical  charge  transfer  between  the  top  and  the  bottom 
surfaces  of  the  slab.  For  the  Al-terminated  (0001)  surface. 


both  the  regular  surface  and  the  surface  model  with  one  N 
atom  adsorbed  at  the  H3  hollow  site  are  considered.  This  N- 
adatom  model  was  shown  previously  to  be  the  most  stable 
structure  under  N-rich  conditions.1161 

The  main  adsorption  configurations  of  atomic  oxygen 
on  Al-terminated  A1N(0001)  surface  are  shown  in  Figure  4, 
panels  a-c.  The  nomenclature  used  to  identify  these 
configurations  is  taken  from  ref.  16.  Among  the  sites  at 
hollow  H3,  hollow  T4  and  atop  sites[161,  the  adsorption  at 
H3  site  is  the  most  stable  with  a  binding  energy  of  204.7 
kcal/mol  with  respect  to  the  isolated  surface  and  an  isolated 
O  atom.  The  adsorption  energy  of  the  T4  site  decreases  to 
168.6  kcal/mol  due  to  the  repulsive  interaction  between  the 
adsorbed  O  atom  and  the  N  atom  in  the  second  surface 
layer.  For  the  most  stable  adsorption  site  H3,  increasing  the 
coverage  to  0.5  ML  decreases  the  adsorption  energy  to  171 
kcal/mol  per  oxygen  atom.  At  a  coverge  of  1.0  ML,  the 
adsorption  energy  drops  further  to  124.3  kcal/mol.  In  the 
case  of  the  N-adatom  model  (Figures  4d  and  4e),  the 
binding  energy  of  the  O  atom  decreases  somewhat  but  it 
remains  strong  with  a  maximum  value  of  108.1  kcal/mol 
corresponding  to  the  H3  site  (Figure  4d). 

In  the  case  of  the  N-terminated  AIN  (0001 )  surface, 
adsorption  on  the  Al-adatom  H3  site,  as  illustrated  in 
Figures  4f  and  4g  and  shown  before1161  to  be  the  most  stable 
site  under  N-rich  conditions,  has  been  investigated.  As 
indicated  in  Table  1,  the  corresponding  adsorption  energies 
of  atomic  oxygen  remain  high  due  primarily  to  A1  atom 
bonding.  These  results  indicate  that  under  N-rich 
conditions,  the  O  atom  can  still  interact  strongly  with  A1 
atoms  on  the  AIN  surface,  independent  of  A1  or  N 
termination  of  the  surface. 

As  an  extension  of  prior  studies  of  the  interactions  of 
nitromethane  (NM)  on  the  a-Al203(0001)  surface,1171  the 
adsorption  properties  of  NM  on  A1N(0001)  surface  have 
been  calculated.  The  results  presented  in  Figure  4h  and  in 
Table  1  indicate  the  existence  of  a  stable  adsorption 
configuration  with  a  binding  energy  of  133  kcal/mol. 
Molecular  bonding  in  this  case  takes  place  by  formation  of 
both  Al-O  and  Al-N  bonds.  From  the  initial  adsorbed 
configuration  (see  Figure  4h)  the  NM  molecule  can  easily 
dissociate  by  overcoming  a  barrier  of  about  8  kcal/mol  with 
formation  of  adsorbed  CH3NO(a)  +  O(a)  fragments.  In  a 
second  reaction  step,  the  remaining  O  atom  of  CH3NO  can 
be  eliminated  with  formation  of  CH3N(a)+20(a).  The 
overall  reaction  is  highly  exothermic  by  161  kcal/mol. 
These  results  clearly  indicate  that  the  A1  terminated  AIN 
surface  can  also  function  to  efficiently  adsorb  and  dissociate 
the  NM  molecule. 
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4.  Conclusions 

A  systematic  study  of  a  series  of  (a)  P-N-C 
heterocycles  containing  azide  ligands,  (b)  salts  consisting  of 
the  [CN7]‘  anion  paired  with  hydrazinium,  ammonium  and 
aminoguanidinium  cations,  and  (c)  compounds  containing 
both  an  oxygen-rich  anion  (dinitramide)  and  a  nitrogen-rich 
cation  (l,5-diamino-4-methyltetrazolium)  interacting  with 
the  aluminum  (111)  surface  have  been  performed  using 
plane-wave  DFT  calculations  in  conjunction  with 
aluminum  surface  slab  models.  These  calculations  indicate 
that,  in  the  majority  of  cases  molecular  adsorption  takes 
place  by  formation  of  multiple  Al-N  or,  in  the  case  of  O- 
containing  compounds,  Al-O  bonds.  For  the  ensemble  of 
molecular  systems  considered  here,  significant  variation  in 
binding  energies  has  been  predicted,  depending  on  the 
number  of  surface  bonds  involved  as  well  as  the  surface  site 
and  the  molecular  orientation.  Additionally,  adsorption  is 
found  to  lead  to  significant  molecular  deformations  and 
upward  shifts  of  some  of  the  surface  aluminum  atoms.  For 
the  high-nitrogen  species  investigated,  dissociative 
adsorption  readily  occurs  leading  to  desorption  of  N2 
molecules.  Atomic  oxygen  binds  strongly  to  the  AIN 
surface.  In  such  cases  the  presence  of  surface  A1  atoms  or 
Al-adatoms  is  sufficient  to  lead  to  large  adsorption  energies. 
Dissociation  of  NM  was  found  to  be  highly  exoenergetic 
and  to  take  place  with  a  small  activation  energy  on  the  Al- 
terminated  (0001)  AIN  surface. 
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Table  I.  Calculated  adsorbate  binding  energies,  in  kcal/mol,  for  various  systems  analyzed  in  this  study  (indicated  in 
parentheses)  together  with  the  corresponding  supercell  model  used  in  simulations. 


Adsorbate 

Binding  Energies  (kcal/mol) 

Supercell  Model 

c2n15p 

60.1  (la),  115.8  (lb),  138.7  (lc),  152.4  (Id) 

6x6 

cn18p2 

55.8  (le),  139.9  (If),  157.8  (lg),  180.4  (lh) 

6x6 

N21P3 

49.5  (li),  80.9  (lj),  110.3  (lk),  149.3  (11) 

6x6 

[N2H5][CN7] 

41.0  (2a),  63.0  (2b),  65.5  (2c),  80.9  (2d) 

4x4 

[NH4][CN7] 

14.0  (2e),  40.9  (2f),  44.4  (2g),  93.5  (2h) 

4x4 

[CN4H7][CN7] 

24.8  (2i),  34.7  (2j),  54.5  (2k),  91.6  (21) 

4x4 

DMTD 

20.4  (3a),  121.9  (3b),  127.5  (3c),  242.9  (3d) 

4x4 

0/A1N(0001) 

204.7  (4a),  168.6  (4b),  150.4  (4c) 

2x2 

O/N  adatom  H3/A1N(0001) 

108.1  (4d),  98.7  (4e) 

2x2 

O/Al  adatom  H3/A1N  (OOOl ) 

100.5  (4f),  80.0  (4g) 

2x2 

NM/A1N(0001) 

133.0  (4h) 

3x3 

This  page  is  Distribution  A:  approved  for  public  release;  distribution  unlimited. 


Figure  2.  Side  views  of  representative  binding  configurations  of  [N2H5][CN7]  (a-d),  [NH4][CN7]  (e-h)  and  [CN4H7][CN7]  (i-1). 
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Figure  3.  Pictorial  views  of  the  adsorption  configurations  of  DMTD  on  Al(lll)  surface. 
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Figure  4.  Side  and  top  views  of  the  oxygen  adsorption  configurations  on  A1N(0001)  (a-c),  N-adatom  A1N(0001)  (d-e)  and  Al- 
adatom  AIN  (000 1)  (f-g).  Configurations  (a-c)  correspond  to  adsorption  at  hollow  FI3,  hollow  T4  and  on -top  sites,  on  A1 
terminated  AIN(0001)  surface.  Panels  d)  and  e)  correspond  to  O  adsorption  at  H3  and  T4  sites  in  the  presence  of  a  coadsorbed  N 
atom  at  H3  site.  Panels  f)  and  g)  correspond  to  O  adsorption  on  N  terminated  Al(000-1)  surface  when  a  coadsorbed  A1  is 
presented  at  a  FI3  hollow  site.  Panel  h)  indicates  the  minimum  energy  potential  pathways  for  dissociation  of  nitromethane  on 
A1N(0001)  surface  leading  to  formation  of  NCH3(a)  radical  and  two  adsorbed  O  atoms  . 


This  page  is  Distribution  A:  approved  for  public  release;  distribution  unlimited. 


